Cook Inlet, Alaska 



R. L. Bruhn* 
W.T. Parry 
M. P. Bunds 



Department of Geology and Geophysics, University of Utah, Salt Lake City, Utah 841 12, USA 



ABSTRACT 

Faulting controls fluid migration within 
transpressive fault-propagation folds in the 
Cook Inlet forearc basin of south -central 
Alaska. Na-Ca-G brine migrates put of Meso- 
zoic rocks through reverse and oblique-slip 
faults into the cores of anticlines, where the 
fluid spreads laterally outward into lower 
Tertiary strata By flow through cross faults 
and peiineabie beds. Precipitation of zeolite 
and carbonate cemerif a^ ■ 
permeability of folded Redding and faults. 
Zeolite minerals are formed by chemical re- 
actions between Na-Ca^CT brine and -j£e«ii-'- 
mentary rocks. Carbonate minerals are pre-- 
cipitated when Na-HC0 3 connate fluid in 
the Tertiary section reacts with rocks during 
diagenesis, and ; by mixing of migrated 
Na-Ca-Cl brine with the Na-HC0 3 pore 
fluid. Carbonate cement is also precipitated 
by fluctuations in PC0 3 during faulting and 
jointing. ' ; t 

High fluid pressure is encountered while 
drilling through lower Tertiary and Meso- 
zoic strata in some anticlines. Hi gh-pressure 
fluid is contained within porous beds that 
are intercalated with strata cemented by 
carbonate and zeolite minerals. Zeolite and 
carbonate cemented beds retard the dissipa- 
tion of high fluid pressure, and channel fluid 
flow parallel to bedding within the anti- 
clines. High fluid pressure may be generated 
by several processes, acting either alone or 
together. The evidence for fault-controlled 
migration of fluid out of the basement sug- 
gests that volumetric strain related to defor- 
mation is most important, but may be aug- 
mented by dynamo-thermal metamorphism, 
sedimentary compaction, alteration of or- 
ganic-rich rock and hydrocarbons, and pos- 
sibly glacial loading. 



Keywords: fluid pressure, forearc basins, 
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Fluid pressures significantly above or below 
hydrostatic occur in a variety of tectonic settings, 
from stable regions to actively aeforming plate 
margins (Ortoleya et aL, 1995; Powjey, 1990; 
Hunt, 1 990). Abnormal pressure is found in com- 
partments funded by low-permeability rocks 
that inhibit pressure equilibration between adja- 
». ; cent regions. Clay-rich strata are. most often cited 
as potential seals (Denting, 1994)^ but recent 
studies have proposed that seals may also be 
formed by diagenetic and hydrbthermal minerals 
(Ortoleva et aL, 1995). 

The goals of this study are to detennine the 
distribution of fluid pressure, to determine the 
origin of cement and vein minerals that inhibit 
fluid flow, and to evaluate mechanisms for fluid . 
.{ . migration and fluid-pressure generation within 
• the upper, or northeastern, part of Cook Inlet 
basin (Fig. 1). Active tectonics; together with ex- 
tensive subsurface data collected during more 
than 30 years of hydrocarbon exploration and 
production, make the basin an attractive natural 
laboratory for studying the interplay between 
geodynamics, fluid-pressure generation, fluid 
migration, and diagenesis. Some previous studies 
conclude that high-pressure fluid is trapped by 
calcite-cemented seals, which cut across strati- 
graphic boundaries in oil field anticlines and 
occupy steeply dipping faults (Powley, 1990; 
Hunt, 1990). Alternatively, others conclude that 
there is little evidence for anomalous fluid pres- 
sure in seal-bounded compartments (Franks and 
He-Zhiyong, 1995); A resolution to this problem 
is important for several reasons. Migration of 
hydrocarbons is partly controlled by the spatial 
and temporal distribution of fluid pressure, but 
information on the fluid pressure regime has not 
been integrated into published petroleum system 



models (Kirschner and Lyon, 1973; Magoc 
1994). Faults in Cook Inlet basin may genen 
damaging earthquakes in the most populated ai 
industrialized part of Alaska (Detterman et a 
1974; Haeussler, in press). Studying spatial vai 
ations in fluid pressure may help to identify mi 
chanically unstable faults, with high potential ft 
earthquake generation. If caltite seals' cut acros 
stratigraphy and structure (Powley, 1990; Hun' 
1990), then fault dccollernent may be controlle 
by the thermodynamics of mineral dissqlutioi 
and precipitation rather than by rhechanica 
anisotropy imparted by sedimentary layering 

f RESEAKCT METHODS t . V 

.. Research work included making observations 
. and collecting samples from outcrops and drill 
cores; compiling and analyzing well log, fluid 
pressure, and rock permeability measurements; 
analyzing geochemistry of cement and vein min- 
verals; and geochemical modeling of cfjerriical 
reactions among basin fluids and rocks ARCO 
Alaska Inc. provided access to seismic reflection 
data that were used to construct geologic cross 
sections of selected anticlines. Construction of 
the cross sections and interpretations relevant to 
earthquake generation were discussed in detail 
by Haeussler et aL (in press). In this study we fo- 
cus on the relationship between faulting, folding, 
and fluid migration. 

TECTONICS OF UPPER 
COOK INLET BASIN 



Cook Inlet is located in a forearc basin bounded 
by the Chugach and Kenai Mountains to the 
southeast and the Alaska Range and Aleutian vol- 
canic arc to the northwest (Fig. 1 ; Dickinson and 
Seeley, 1979). The geodynamics of Cook Inlet 
basin are complex because of the protracted his- 
tory of deformation associated with subduction 
and rnicroplate collision (Plafker et aL, 1994). Re- 
gional transpression is caused by mechanical cou- 
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pling between the North American and Pacific 
plates, and collision of the Yakutat block at the 
eastern end of the subduction zone. The depth to 
the subducting Pacific plate is 50-60 km beneath 
the basin's center and deeper on its northwest side 
(Page etal., 1991). 

Large faults bound three of the basin's mar- 
gins: the Casde Mountain fault at the northeast 
end of the basin, the Bruin Bay fault in the 
foothills of the Alaska Range, and the Border 
Ranges fault along the northern flanks of the 
Chugach and Kenai Mountains (Fig. 2; Grantz, 
1966; Pavlis, 1982; Plafker et al. a 1994)i Paleo- 
zoic and Mesozoic rocks in the Alaska Range are 
part of an accreted arc and microcontinent that 
collided with North America during Cretaceous 
time (Silberling et al., 1994). These rocks extend 
southeastward beneath Cook Inlet basin and out- 
crop as fault-bounded slivers and klippen along 
the flanks of the Chugach and Kenai Mountains. 
Mesozoic sedimentary rocks in the Chugach and 
Kenai Mountains are part of a vast accreted com- 
plex that was faulted against and thrust beneath 
the southern edge of the continental margin. Ter- 
tiary strata in Cook Inlet t>a$irV are therefore de- 
posited on a polyglot basement of Paleozoic and 
Mesozoic metamorphic and igneous rocks, and 
Jurassic through Cretaceous marine strata. 

Cook Inlet basin formed as Paleocene streams 
deposited strata over deformed Mesozoic and 
older rocks (Fig. 3; Calderwood and Fackler, 
1972). Lower Tertiary strata and volcanic rocks 
were deformed by regional strike-slip faulting 
and folding during Eocene and early Oligocene 
time (Barnes and Payne;* i 956; Clardy, 1974; 
Starnataikbs et aL t 1988* 1989; Little, 1991). 
Transprcssion formed a complex structural terrain 
of folds, faults, eroded horst blocks, and adjacent 
grabens filled with synorogenic fluvial deposits. 
Conglomerate and coarse sandstone were de- 
posited over much of the basin as deformation 
waned, forming the Hemlock, Bell Island, and 
Tsadaka Formations (Calderwood and Fackler, 
1972; Clardy, 1974). Deposition of lower Mio- 
cene fluvial strata buried this middle Tertiary 
structural terrain prior to the onset of late Mio- 
cene to Holocene deformation (Kirschner and 
Lyon, 1973; Calderwood and Fackler, 1972). 
Neogene deformation continued along parts of 
the Bruin Bay and Castle Mountains faults, and 
initiated fault-propagation folding throughout the 
basin (Fig. 4; Kirschner and Lyon, 1973; Boss 
et al., 1976). The anticlines contain prolific oil 
and gas reservoirs (Kirschner and Lyon, 1973; 
Magoon, 1994). 

The anticlines are cored by reverse to oblique- 
reverse slip faults that root in the Mesozoic and 
older basement rocks and cut upward through 
lower and middle Tertiary strata (Fig. 4; Kirschner 
and Lyon, 1973; Haeussler et al., in press). The an- 
ticlines are also cut by numerous cross faults that 




Figure 1. Tectonic map of the southern Alaskan plate margin showing the location of ( 
Inlet Basin. The cross section shows the distribution of earthquakes in the BeniofT zone 
North American plate in the Cook Inlet region for a 10 yr period. Earthquakes are proje 
onto a west-northwest-east-southeast vertical plane. Earthquake data provided by John L 
U.S. Geological Survey seismicity catalog. 



offset fold hinges, and subdivide oil and gas fields 
into distinct hydrologic units (Alaska Oil and Gas 
Conservation Commission, 1 994). Kirschner and 
Lyon (1973) proposed that displacement on cross 
faults was mostly normal slip, but both outcrop ex- 
posures and structural contours of cross-faulted 
folds indicate predominantly strike-slip displace- 
ment (Fig. 4A; Barnes and Payne, 1956; Bruhn 



and Pavlis, 1981; Alaska Oil and Gas Conse; 
tion Commission, 1994). 

Carbonate and zeolite veins and cement oc 
in both large reverse-slip faults and smaller cr 
faults where observed in lower Tertiary and ol 
strata (Fig. 5). The veins formed during multi 
episodes of fluid flow and deformation. Evide: 
for mineral precipitation during faulting indue 
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Selected Oil & Gas Fields 
BR: Beluga River 
NC: North Cook Inlet 
GP: Granite Point 
TB: Trading Bay 
MGS: Middle Ground Shoal 
MC: McArthiir River 
SR: Swanson River 
BC: Beaver Creek 
KC: Kenai - Cannery Loop 



Inset 2 



SCU-33-33 



Beaver Creek 
Fold JT 



10 km : 



Figure Z Tectonic map of upper Cook Inlet showing large faults, prominent anticlines and oil 
and gas fields. Insets 1 and 2 enclose North Cook Inlet and Swanson River anticlines. Numbers 
indicate locations of wells listed in Table 1. 



multiple sets of crosscutting veins, and mineral 
fibers elongated parallel to fault striae in both 
drill core and outcrop. Faults in late Miocene and 
younger rocks are discrete zones of cataclasite 
and gouge with few vein minerals. 

The distribution of fluid pressure within Cook 
Inlet basin must be considered in context of the 
history of hydrocarbon generation, migration, and 
trapping. Reservoirs in the upper Miocene Beluga 
Formation and Pliocene Sterling Formation pro- 
duce gas generated by biogenic alteration of coal 
beds (Claypool etal., 1 980). Oil and gas produced 
from the lower Miocene Tyonek Formation and 
older strata are primarily derived from thermal 
maturation of organic-rich, marine source beds 
Ln the Middle Jurassic Tuxedni Group (Magoon 
and Claypool, 1981). Hydrocarbon migration oc- 



curred in two phases, according to Magoon and 
Claypool (1981) and Magoon (1994); Hydrocar- 
bons first migrated into traps along the angular 
unconformity between Mesozoic and lower Ter- 
tiary rocks. These traps were breached by defor- 
mation and erosion during middle Tertiary time. 
Deposition of Miocene strata was followed by re- 
newed migration during late Miocene to Holo- 
cene faulting and folding Hydrocarbons became 
trapped in the crests of growing anticlines during 
this latest phase of migration. 

LITHOLOGY OF TERTIARY AND UPPER 
MESOZOIC SEDIMENTARY ROCKS 

Tertiary strata are sandstone, conglomerate, 
siltstone, and coal deposited in meandering and 



Calderwood and Fackler, 1972). These dep 
are as much as 9 km thick in the deepest pa 
Cook Inlet basin, but Miocene and youi 
strata thin over the crests of anticlines bee; 
of erosion and stratal offlap caused by fok 
(Kirschner and Lyon, 1973; Boss et al., IS 
Haeussler et al., in press). Detrital grains are 
termediate composition volcanic, metamorpi 
and plutonic lithic fragments, and quartz, j 
gioclase, and less common K-feldspar, deri< 
from erosion of the surrounding mountai 
Middle to upper Mesozoic strata are mos 
lithicwacke, siltstone, mudstone, shale, a 
minor limestone deposited in a marine envirc 
ment over a Jurassic arc and older basemei 
However, fluvial strata and coal beds also occ 
in Upper Cretaceous deposits (Magoon et a 
1980). 

Cement and vein minerals decrease the porosr 
and permeability of sedimentary rocks and fault 
forming aquitards to fluid flow. Cement minera 
include clay, quartz, chlorite, calcite, sideritt 
dolomite, and zeolites (Hayes,. 1979; Lyle an 
Morehouse, 1977; Fisher and Magoon, 1978; 
Secondary carbonate and zeolite cements ar 
abundant in lower Miocene and older sedimen tar 
rocks, where calcite and siderite fill pores and re 
place original feldspar grains, clay, chlorite, anc 
quartz. Laumonute and prehnite, as well as car- 
bonate minerals, seal fractures and cement pores in 
the Paleocene West Foreland Formation and oldei 
rocks throughout much of the basin. 

The carbonate- arid zeolite-cemented rocks 
form low-permeability, lenticular stringers and 
tabular beds in tercalated between more porous and 
permeable beds from which cement and frame- 
work grains were partly dissolved. The lateral di- 
mensions of secondary * cemented and secondary 
porosity-enhanced beds could not be established 
between wells. However, siderite and calcite- 
cemented sandstone and conglomerate beds, me- 
ters to tens of meters thick, extend for tens to hun- 
dreds of meters along strike in the mining high 
wall of Paleocene coal measures in the Wishbone 
Hill fold in the Matanuska Valley (Fig. 5). This ex- 
humed fold is structurally analogous to folds in the 
subsurface of upper Cook Inlet basin. 

FLUID-PRESSURE MEASUREMENTS 

Methods 

Information on fluid pressure was compiled 
from published summaries of initial reservoir pro- 
duction pressure (Alaska Oil and Gas Conserva- 
tion Cornmission, 1994), from drill stem test 
(DST) and repeat formation test (RFT) records, 
and by interpretation of well log petrophysics 
measurements obtained from reports and files of 
the Alaska Oil and Gas Conservation Commis- 
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. Figure 3. Stratigraphy with annotated tectonic history of upper Co 
basin (upper column) and the adjacent Matanuska Valley (lower c 
(see Fig. 2). 



sion. We selected pressure records from explo- 
ration wells after determining that reported pres- 
sures were not in error because of leaking packers 
or other problems during testing (Table 1); 

Initial Reservoir Pressure. Initial fluid pres- 
sures in Cook Inlet gas and oil reservoirs were 
tabulated and published by the Alaska Oil and 
Gas Conservation Cornmission (1994; Fig. 6A). 
Initial pressures generally increase systemati- 
cally with reservoir depth at a rate of -10-1 1 
MPa/km (Fig. 6 A). Initial pressure is between 
10% and 20% above hydrostatic in several shal- 
low gas reservoirs in the Beluga, Sterling, and 
l^onek Formations, and in oil reservoirs in the 
Hemlock Conglomerate in the Swanson River 
oil field. Few production reservoirs are located 
below a depth of 3 km, the proposed transition 
depth from normal to high fluid pressure in the 
basin proposed by Hunt (1990). Evidence for 
fluid pressure at greater depth must be gleaned 
from fluid-pressure measurements and petro- 
physics logs collected in exploration wells. 

Fluid-Pressure Measurements, DST and 
RFT results from exploration wells are listed in 
Table 1 and plotted with respect to depth in Figure 
6B. Several aspects of these fluid pressure mea- 
surements are noteworthy. 

1. Fluid pressure is mostly hydrostatic above 
3 km depth, but becomes greater than hydro- 



static in some deeper intervals, as discussed by 
Hunt (1990). ■:■ 

2. Fluid pressure fluctuates with depth rather 
than increasing monotpnically along either 
hydrostatic or abnormally high pressure gradi- 
ents. Consider the Sunfish #3 exploration well in 
the North Cook Inlet anticline (Fig. 2; Table 1), 
for example. Fluid pressure in the Tyonek For- 
mation is hydrostatic until the C-sands are en- 
countered at 3400 m, where fluid pressure is 
50 Mpa, or 33% greater than hydrostatic (Fig. 6B; 
Table 1). Fluid pressure beneath the C-sands de- 
creases to 43 Mpa, or nearly hydrostatic pressure 
at 3750 m (Table 1). 

3. Fluid pressure varies between wells in the 
same stratigraphic interval. For example, fluid 
pressure in the lower Tyonek Formation C-sands 
is significantly different in the Sunfish #1 and #3 
wells, which both penetrate the North Cook Inlet 
anticline (Eg. 2; Table 1). Fluid pressure in the 
C-sands is approximately hydrostatic in the Sun- 
fish #1 well, but is 33% above hydrostatic in the 
Sunfish #3 well located -10 km to the north in 
the same anticline. 

4. Fluid pressure is greater than hydrostatic in 
Mesozoic rocks in the Swanson River anticline. 
High fluid pressure was measured in the Upper 
Cretaceous Matanuska Formation and the Upper 
Jurassic Naknek Formation (Table 1). 



Petrophysics Measurements. Sonic > 
logs provide additional evidence for the ^ 
distribution of overpressure in the core 
Swanson River anticline/ where DSTs ir 
high fluid pressure in Mesozoic rocks (T 
wells SCU-33-33 and SCU-22-32). The ii 
transit time (i l l ) of compressiona] wave* 
ing through rocks next to the well bore u 
decreases with depth because porosity is re 
by compaction, causing sonic velocity to iro 
(Serra, 1986). However, rocks with high 
pressure are undercompacted. Undercomp. 
results in lower compressional wave veloci 
increased ITT. This is illustrated by the 
velocity log in the SCU-33-33 well (Fij 
where we have plotted ITT in sandstone 1 
The ITT decreases systematically along a roi 
linear trend on the semilogarithmic plot, 
depth of -—4000 m. Below -4000 m the ITT 
tuates about a nearly vertical trend in iin 
compacted sandstone beds that are intercai 
with mudstone and shale. We interpret this : 
velocity anomaly to be caused by overpressi 
Mesozoic sandstone in the core of the fold. 

PERMEABILITY MEASUREMENTS 

Permeability measurements are compi 
from well test and commercial laboratory rcpi 
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Figure 4. (A) Structural configuration of Hemlock conglomerate in the Swanson River anti- 
cline constructed from structural contour maps prepared by oil field operators and exploration 
companies. Lateral and vertical offsets of the fold crest are caused by cross faulting, perhaps like 
that in the Wishbone Hill fold (Fig. 5). (B) Geologic cross section of the North Cook Inlet anti- 
cline showing reverse or oblique-reverse faults in the core of the anticline. 



on file with the Alaska Oil and Gas Conservation 
Commission. The permeability measurements 
are mostly from sandstone and siltstone in poten- 
tial reservoir strata, rather than lower permeabil- 
ity muds tone and shale. The equipment used by 
commercial laboratories also restricts measure- 
ments of low-permeability rocks. The lower per- 
meability limit for most commercial laboratory 
tests was ~10~ 17 m 2 , with the exception of core 
samples from the Sun Fish #3 well, where the 
instrument detection limit was -10~ 18 m 2 . Only a 
few tests are made at elevated confining pressure; 
most are air-permeability tests at atmospheric 
pressure. Here we discuss core plug measure- 



ments from the Sun Fish #3 well in the North 
Cook Inlet anticline. Several Tertiary and Meso- 
zoic formations were tested in this well; some of 
the tests were done at elevated confining pres- 
sure, and in several intervals core plugs were 
drilled both parallel and normal to bedding. 

Permeability is reported for strata in the Meso- 
zoic Naknek and Matanuska Formations, and the 
overlying West Foreland, Hemlock Conglomer- 
ate, and Tyonek Formations in the Sunfish #3 
well. Permeability (k) spans several orders of 
magnitude, ranging from the lower detection limit 
ofk=* lCT I8 m 2 to a high of k = 10" ,3 m 2 in the 
most permeable samples. Permeability measured 



greater than perpendicular to bedding. Co 
grained rocks are more permeable than 
grained rocks, with the exception of intervals 
secondary carbonate and zeolite cement In I 
latter intervals, the permeability of conglom* 
and sandstone decreases markedly, and in s 
cores conglomerates are almost as impermes 
or even more so, than adjacent siltstone, crea 
alternating bands of higher and lower permez 
ity that mimics vertical variations in bedding. ' 
vertical distribution of secondary cement flu< 
ates markedly between adjacent beds, and e* 
within the same bed. This variability is detec 
either by visual inspection of the core, or by n 
: ; ing that the average grain density varies invers. 
with permeability because of high-density calc 
and siderite cement (Fig. 8). 

The permeability of packed off intervals se 
era! meters to tens' of meters long was determin 
as part of DSTs and RFTs in several wells. T 
permeability ranges from k < 10" 18 m 2 in carbo 
ate cemented intervals to much greater perme 
bility in less, cemented sandstone and congli 
merate (k « HTJ 4 m 2 -tp iO" 12 m*> For exampl 
k £ lO" 18 m 2 in lower TVonekFonnariob sane 
stone and siltstone in the Big Lake #1 explcratio 
well southeast of the Castle Mountain fault zone 
even though the formation has been uplifted sev 
eral kilometers by folding and faulting. 

FLUID COMPOSITION IN MESOZOIC 
AND TERTIARY ROCKS 

.... ;, ? > ; .... 

Two chemically dis tinct pore fluids are preseni 
in Cook Inlet basin below a shallow* fresh- to 
brackish-water transition (Table 2; McGee, 1977; 
Franks and He-Zhiyang, 1995). Oil and wet gas 
in lower Tertiary strata are associated with a 
Na-Ca-CI brine containing 17 000 ppm to 35000 
ppm dissolved solids and relatively heavy 5 18 0 
and 5D. This fluid is altered seawater derived 
from Mesozoic marine source beds during dia- 
genesis, and produces high-conductivity anom- 
alies around the crests of some anticlines (McGee, 
1977). Dry gas in late Miocene and younger 
reservoir rock is associated with an Na-HC0 3 
brine with 3500-17000 ppm, but typically be- 
tween 3500 and 5000 ppm, total dissolved solids 
and lighter 5 ,8 0 and 5D than the Na-Ca-CI brine. 
This fluid is derived by alteration of fresh water 
trapped and buried in the Tertiary strata (Franks 
and He-Zhiyang, 1995). 



STABLE ISOTOPE COMPOSITION OF 
CALCITE CEMENT AND VEINS 

Carbon and oxygen isotopes were measured 
in calcite cement and veins in samples of Meso- 
zoic and Tertiary rocks collected from drill 
cores. Veins were also collected from strike-slip 
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figure 5. Geologic map of the Wishbone Hill fold, Matanuska Valley (Fig. 2). Resistant ridge of calcite and siderite veins and cement o 
the Eska fault, a left-lateral cross-fault that offsets Paleocene and Eocene strata in the Wishbone Hill fold. Qac — Quaternary glacial anc 
deposits; Qls — Quaternary landslide; Tar — Arkose Ridge Formation (Paleocene); Tc — Chickaloon Formation (Paleocene-Eocene); TV- 
bone HuT Formation (Eocene); Tt— Tsadaka Formation (Oligocerie); and Km — Matanuska Formation (Upper Cretaceous). 



faults in the Wishbone Hill fold and Castle 
Mountain fault zone in the Matanuska Valley 
(Figs. 2 and 5). The samples show varying de- 
pletion of 5 13 C from about 0%o to -25%o relative 
to Peedee belemnite (Fig. 9). Depletion of 8 ,3 C 
below -5%o is presumably caused by the influ- 
ence of an organic carbon component. 

There are three distinct trends in the 5 l8 0 data 
(Fig. 9). Cements from Tertiary rocks are rela- 
tively light, cements from Mesozoic rocks are 
relatively heavy, and veins from faults are inter- 
mediate in oxygen isotopic composition. The 
5 ls O values of calcite cements presumably 
reflect isotopic equilibration between connate 
waters present in the rocks at about 80 °C. Frac- 
tionation between calcite and water at 80 °C is 
\9.5%o; thus the +20%o values for Mesozoic 
cements represent equilibration between calcite 
and water with 5 ls O = 0%o, and +6%o calcite in 
Tertiary rocks equilibrated with -13.5%o water. 
Several calcite cement samples from lower Ter- 
tiary strata have intermediate and heavy S I8 0, 
presumably due to either brine mixing or equili- 



bration of the rock with fluid that migrated up- 
ward from the Mesozoic basement. The isotopic 
compositions of these inferred fluids match those 
of the Na-Ca-CI and Na-HCO, fluids sampled in 
oil wells (Table 2). Intermediate 5 ls O values in 
cement and veins could result from either mixing 
of these two fluids, or by fluid-rock reactions. 
The efficacy of these processes was investigated 
by geochemical modeling. 

GEO CHEMICAL MODELING OF 
ZEOLITE AND CARBONATE 
PRECIPITATION 

Precipitation of secondary calcite and zeolite 
cement and vein minerals reduced the permeabil- 
ity of lower Tertiary and older strata, and reverse- 
and strike-slip faults. Zeolite and carbonate 
cement may have precipitated during mixing of 
Na-Ca-CI and Na-HC0 3 fluids, by reactions 
between these fluids and minerals in the rocks, 
and by changes in the partial pressure of C0 2 in 
pore fluid. These processes were modeled with 



the speciation and reaction path pre 
SOLVEQ and CHILLER (Reed, 1 982). Te 
ture was fixed at 80 °C in all simulations, t 
lent to the temperature at a depth betweei 
and 4 km in upper Cook Inlet. Starting coi 
tions of the two fluids were chosen witr 
measured range of Na-Ca-CI and Na- 
brines (Franks and He-Zhiyong, 1995; Tal 
Mixing Na-Ca-CI brine with Na-HCO, 
precipitated both calcite and siderite (Fig. 
Example results are cited for reference: mixii 
kg of Na-Ca-CI brine with 1.0 kg of Na-f 
brine produced 1 .6 g of calcite or 0.59 cm 3 c 
cite per 1.5 kg of water. A total of 2.5 x 10 5 
water, or 255 m 3 , is required to fill the pore vo 
in 1 m 3 of rock with an initial porosity of 

Chemical reactions of Na-Ca-CI and Na-H 
fluids with a sandstone consisting of 50% qu 
40% plagioclase (An 42 ), 5% K-feldspar, 
chlorite, and 2% muscovite were also simulc 
The simulations were carried to albite saturat 
although plagioclase continued to react. Mu- 
vue, quart2, calcite, prehnite, laumontite, 
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Well Name 


Reference 


(») 


ueoin 








Sunfish 1 (SF1) 


Inset 1 


9562 
10423 


2898 
3158 


4295 
4778 


30 
33 


Tyonek 
Tyonek 






10856 


3290 


5188 


36 


Tyonek 






10860 


3291 


5208 


36 


Tyonek 






•?707 

O/Ui 




43 


Tyonek 


Sunfish 3 (SF3) 


Inset 1 


11253 


3410 

O / CO 


7200 
6276 


50 

43 ' 


Tyonek 
Tyonek 


N.Cook Inlet St (NCI1) 


Inset 1 


4033 
11015 


1222 


1691 
7579 


12 
52 


Sterling 
Tyonek 


SRS Slate 1 (SRS1) 


Inset 1 


15040 
16050 


4558 
4ao4 


12893 

1 1d OA 


89 
93 


TVonek 
W. Foreland 


Big Lake Test 1 


1 


4553 

CCOQ 
DDOO 


1380 

ID/O 


2602 
3187 


18 

22 


Bell Island 
Tyonek 


Naptowne Unit 24-8 


6 


6064 
6333 
6444 


1838 
1919 


2545 
2648 


18 
18 


Sterling 
Sterling 






1953 


2738 


19 


Beluga 






9421 


2855 


4775 


33 


Tyonek 






9445 






33 


Tyonek 


SCU-22-32 


insei d. 




4273 


9120 


63 


Naknek 


SCU-33-33 


inset £ 




4091 


8000 


55 


Matanuska 


Wolf Lake 1 (WL1) 


msec « 




4052 


8000 


55 


Naknek 


Beaver Cr. 5RD (BC5RD) 


insei * 


14953 


4531 


9250 


64 


Tyonek 


BeaverCr.2(BC2) 


Inset 2 


14900 


4515 


7135 


AO 


nollUOCK 




15130 


4oao * 


6959 


48 


Hemlock 


Beaver Cr. 8 (BCB) 


Inset 2 


9233 
9248 


2798 : 

:..2o02 


4600 

AAIO • 


32 
31 


Tyonek 
Tyonek 


Granite Point #1 


3 


8665 
9230 


'.•'2626' : ''- 
2797 


' 4304' . 
4348 


30 
30 


Tyonek 
Tyonek . 






11000 




49 


Hemlock (?) 


MGS 17595 




5540, 
5942 \ 


1679 
180i 


26B3 
2687 


18 
19 


N.A. 
N.A. 






6053 


1834 


2775 


19 


N.A. 






6079 


184? 


2803 


19 


N.A. 






7105 


2153 


3735 


26 


N-A. 






7700 


2333 


3785 


26 


NJV. 






8100 


2455 


3870 


27 


NA 






8400 , 


2545 


3745 


26 , 


NA. . . 






8633" 


2616 


4286 


30' 


N.A. " 






8835-M 


2677 




• . 28 


HA. rjtif-v 


MGS 18745 


4 


1Q4O0 


• 01 CO. 


;*„' 4268 


• 29 ' • 


NA: ^ " v 


TVonek St. 18742 


*> 
c 


8270 


2506* 


■ ' 4025^ 


28 




W. Foreland f 


M A' ; 


Q<Ma 
auoo . 


2829' : * 


4475 


'''31'"'*' 


Tyonek ' 


Note: N.A.^-fiot available. Map refen 


3nce refers 


to well location on either main map or inset maps in Figure i . 



chlorite were precipitated by reaction of the 
Na-Ca-Cl brine with sandstone (Fig. 10B). The 
fluid reached saturation with calcite after reaction 
with 0. 3 g of roc); with chlorite after reaction with 
0.8 g of rock, with prehnite after reaction with 1 g 
of rock, and with laumontite after reaction with 
2.9 g of rock. Reaction with 4 g of rock precipi- 
tated 0.05 g of calcite. One m 3 of rock with 10% 
porosity required 120 m 3 of water to completely 
fill pore space with reaction products. 

Quartz, muscovite, calcite, chlorite, and kaolin- 
ite were precipitated in simulated reactions of rock 
with Na-HC0 3 water (Fig. IOC). Reaction of wa- 
ter with 4 g of rock precipitated 0.25 g of calcite, 
and 1 rn 3 of rock required 136 m 3 of Na-HC0 3 
water to fill 10% porosity with reaction products. 

Calcite precipitation was also simulated by de- 
creasing the partial pressure of C0 2 , a process 
that presumably occurs during fracturing. De- 
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creasing C0 2 pressure from 10 bar to 1 bar pre- 
cipitated 0.025 cm 3 of calcite per 1 kg of water. 
Complete sealing of 1 m 3 of rock with 10% 
porosity required precipitation of calcite from 
4000 m 3 of water. Clearly, fluid mixing and fluid 
reaction with rock are more effective at sealing 
porosity than C0 2 loss. 

GENERATION OF ABNORMAL 
FLUID PRESSURE 

The temporal scale for decay of abnormal fluid 
pressure in a seal-bounded reservoir is a function 
of both the permeability and the thickness of the 
seal (Deming, 1994). Fluid pressure begins to de- 
cay in 1 0 5 yr or less if the sealing interval is £1 km 
thick with permeability k > 10" 18 m 2 (e.g., Dem- 
ing, 1994; Neuzil, 1995). This implies that mecha- 
nisms that generate high fluid pressure are either 
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anility or Mesozoic ana lemary tocks in up 
Cook Inlet basin. We group candidate mechanis 
into several categories for discussion: (1) volurn 
ric strain, (2) hydrocarbon generation and alt 
ation, and (3) mineral dehydration. Several ad 
tional processes are eliminated from considerati 
because of the environment in Cook Inlet bas: 
Aquathermal pressuring is unlikely given the ge 
thermal gradient of 20-25 °C /km (Osborne ai 
Swarbrick, 1997). Petrographic observations inc. 
cate little porosity and permeability reduction t 
pressure solution. Topography-driven flow froi 
mountains surrounding the basin was numerical] 
modeled and rejected as a source of abnormal flui 
pressure by He-Zhiyong and Franks (1995). Q 
motic processes are also unlikely to generate nig 
fluid pressure because shale beds are discontinu 
ous in the Tertiary section. 

The one-dimensional hydrodynamic disequi 
librium equation proposed by Neuzil (1995) i 
used to estimate the efficacy of several altemativi 
processes for generating abnormal fluid pressure 

t^kjl a; 

r is hydrodynamic or geologic forcing in units ot 
sr 1 , K is the hydraulic conductivity of either the 
reservoir or the reservoir bounding seal, and L is 
the half-width of the reservoir as measured across 
its smallest dimension. Permeability magnitude k 
- 10~ 7 ATfor brine-filled pore space. High fluid 
pressure may be maintained indefinitely if in- 
equality 1 is satisfied, otherwise reservoir pres- 
sure decays toward ambient pressure in the sur- 
rounding region (Neuzil, 1 995). 

r is either calculated, or discussed qualita- 
tively, depending on available information con- 
cerning a specific rriechanism We conclude that 
..abnormal fluid pressure may be maintained in 
reservoirs several kilometers or more in size if T> 
10- 14 s~\ given that ID" 18 m 2 <£ k <£ 10" 16 m 2 for 
most intervals of zeolite- and carbonate-cemented 
rock in Cook Inlet basin (Table 3). Presumably, 
the permeability of deformed and cemented shale 
and mudstone layers is <10~ !8 m 2 , but these rocks 
were not tested for permeability. Shale and mud- 
stone permeated by zeolite and carbonate cement 
and veins is more common in the Mesozoic sec- 
tion than in the lower Tertiary rocks. 

Process Category 1: Volumetric Strain — 
Sedimentary Compaction, Tectonic Strain, 
and Glacial Loading 

Volumetric strain alters fluid pressure and in- 
duces fluid flow because pore space either con- 
tracts or dilates. 

Sedimentary Compaction. High fluid pres- 
sure may be generated by compaction and thermal 
expansion of pore fluid during subsidence and 
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burial (Bredehoeft and Hanshaw, 1968). Pliocene 
and Quaternary deposits thicken in synclines and 
thin over the crests of anticlines in upper Cook In- 
let basin. Average sedimentation rates are esti- 
mated by sampling the Pliocene (Sterling Forma- 
tion) and Quaternary deposit isopach map of 
Hartman et al. (1972) at 10 km intervals on a 
square grid, and then dividing the sedimentary 
thickness at each grid point by an elapsed time of 
5.3 m.y. Sedimentation rates range from nearly 
zero along the crest of the Middle Ground Shoal 
anticline to =6 x m/yrin the deepest synclinal 
trough, with the average rate dI/d/» 3 x KHm/yr. 
These rates are approximate because the base of 
the Pliocene Sterling Formation is difficult to de- 
tect in some well logs, and the base of the Sterling 
Formation may not mark the Miocene-Pliocene 
temporal boundary at all localities (S. Franks, 
1997, personal commun). 

Hydrodynamic disequilibrium by sedimentary 
compaction is estimated (Neuzil, 1995): 

r=-|[(l + v)/3(l - v)]kV-C t G|dLHr (2) 

where G = geo thermal gradient (20-25 °C/km), 
C m = 1 x 10-* Pa" 1 , C, = lfr 4 C°^, Poisson's ratio 
v = 0.25, y a = 2.3 x 10 4 kg nr 2 s -2 (specific 
weight of fluid saturated reck), and dUdt = sedi- 
mentation rate (m/yr). 

Setting ah upper bound thermal gradient G = 
25 °C/km and dUdi = 6x10"* m/yr in equation 2 
yields T£4 x lO" 15 sH (Table 3). 

Tectonic Strain. Anomalous fluid pressure 
may be generated by volumetric strain during 
faulting and folding (Neuzil, 1995). We consider 
three different tectonic mechanisms: distributed 
deformation, tectonic compaction by shearing in 
fault zones (Blanpied et ah, 1992), and porbelas- 
tic deformation (Scholz, 1990). Consider the case 
where deformation is distributed evenly across 
the basin, with a regional shortening rate of 
-1 cm/yr, or -20% of the plate convergence 
velocity. The linear strain rate ^=-3 x 1CT 15 s" 1 . 
Following Neuzil (1995), we set volumetric 
strain rate Ey = 0.1 and calculate: 
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+ Hemlock Congl. 
w West Foreland Fm. 




-5000 



10 



B 



20 30 40 
Reservoir Press (MPa) 

Pressure (MPa) 
0 50 100 



50 



60 



-1000 - 



/E -2000 - 



3 



-3000 - 



-4000 - 



-5000 L 




Figure 6Y(A) Initial reservoir pressure versus reservoir depth. Solid lines denote hydro: 
pressure gradient for water, estimated hydrostatic gradient for brines in Middle Ground i 
oil and gas field, and a reference line for fluid pressure gradient that is 20% in excess c 
Middle Ground Shoal gradient (MGS). Fluid pressure data are from AOGCC (1994). Mi 
Ground Shoal pressure gradients are from drill stem tests on file with Alaska Oil and Gas < 
seryation Commission. (B) Fluid pressure determined by drill stem and repeat formation 
in Cook Inlet exploration and production wells (see Table 1 for details). Solid line is the rt 
ence hydrostatic gradient Abnormal, high-pressure measurements are enclosed in the poly, 



r=Ka 



(3) 



where £ = 0. 1 is the coefficient relating bulk vol- 
umetric strain to porosity loss (Neuzil, 1995), and 
r « 3 x ICH 7 s~*, two orders of magnitude less 
than forcing by sedimentary compaction. 

Tectonic compaction and consolidation of fine- 
grained fault rocks decreases both porosity and 
permeability, traps fluid, and creates high pore 
pressure (Sibson, 1990; Sleep and Blanpied, 
1 992; Byerlee, 1 993). The high pore pressure may 
trigger rupturing and release high-pressure fluid 
that flows vertically upward and laterally outward 
from the fault zone (Sibson, 1990; Blanpied et aL, 
1992; Byerlee, 1993). 



Consider an active shear zone 10-1 00 m wide, 
and a displacement rate of 1-3 mm/ yr; these are 
appropriate parameters for the Castle Mountain 
fault The upper and lower bound shear strain 
rates for this combination of displacement rates 
and shear-zone widths are: 3.2 x 10~ 12 s"* 1 < y £ 
9.2 x 10~ n s~ ! . If the rate of compaction normal 
to the shear zone is 10% of y, then 

3.2xlO- , V 1 <r<9.2xlO- I3 s- , f 

a volumetric strain rate that is orders of magni- 
tude greater than that estimated for regionally 
distributed shortening strain, and one to two 
orders of magnitude greater than hydrodynamic 
forcing by sedimentary compaction. 

Poroelastic strain creates pressure pulses rang- 
ing from a few tenths to several MPa in com- 



pressional and dilatational lobes around the 
of a ruptured fault (Scholz, 1990). Poroela 
pressure pulses apparently decay within 0. 
100 yr because fluid diffuses into the cour 
rock and equilibrates pressure between the cc 
pressional and dilatational lobes (H. Wa 
1998, personal commun.). This decay time 
much shorter than the estimated recurrer 
interval of >10* yr for reverse and oblique-s 
faults in Cook Inlet anticlines (Haeussler et i 
in press). 

Glacial Loading. High fluid pressure rn, 
be generated in rocks beneath glacial ice shee 
where fluid is trapped in low-permeabilii 
rock, or in reservoirs surrounded by low-pc 
meability rock (Bahr et al., 1994; Thorscr 
1996). T is estimated by substituting the spe 
cific weight of ice (y= 9.0 x 10 3 kg nr 2 s~ 2 ) fo 
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Figure 7. Interval transit time (ITT) in sandstone beds as a function of depth in the SCU-33-33 
weD, in the Swanson River fold; Each gray cross represents the interval transit time in micro- 
seconds per foot of compression*! waves in sandstone adjacent to the well bore. Sandstone layers 
were identified using the spontaneous potential and gamma ray logs following procedures 
described by Serra (1986). Straight line represents the trend of ITT with depth between 1200 m 
and 3200 m found by linear regression. Note the high ITT (low sonic velocity) in a bed just above 
the Tertiary^Mesozoic rock unconformity, and the extensive zone of low velocity below 4000 m. 
The low compressional wave velocity reflects undercompaction that is presumably caused by 
high fluid pressure. t v , . . - : :..£V\ 



that of water-saturated rock and deleting the 
thermal, term in equation 1. Assuming an ice- 
thickehing (accumulation) rate of 1 0~ 2 m/yr to 
10° ra/vr results in I> 1.5 x lO* 1 ^ s" 1 to 1.5 x 
jr 1 , comparable to hydrodynamic forcing 
in fault zones (Table. 3). Cook Inlet basin was 
glaciated during the late Pleistocene as ice 
lobes advanced into the basin from the Alaska 
Range and Chugach-Kenai Mountains, but the 
glaciers retreated from the basin by =«15 k.y. 
ago (Karlstrom, 1964; Schmoll et al. f 1984; 
Reger, 1983). Consequently, any anomalous 
fluid pressure generated by ice loading is relic 
and decaying. 

Process Category 2: Hydrocarbon 
Generation and Alteration 

Conversion of solid kerogen and coal to liquid 
and gaseous hydrocarbons, and thermal alteration 
of oil to gas, may generate abnormal fluid pressure 
(Neuzil, 1995; Barter, 1987; 1990; Osborne and 
S warbrick, 1 997). Middle Jurassic source rocks 
are overmature in the deeper part of Cook Inlet 
basin, however, and the low gas to oil ratio in lower 
Tertiary reservoirs (gas/oil < 1000) suggests that 
gas production is probably not vigorous at present 



(L. Magoon, 1994, and 1998, personal cornrnun.). 
Thermal alteration of coal beds at temperatures be- 
. tween 70 and 120 °C may liberate methane and 
. C0 2 and generate high fluid pressure (Hunt, 
1 996). This process could be occurring in Upper 
Cretaceous and lower Tertiary strata in parts of the 
basin. Notably, alteration of organic material may 
also raise the partial pressure of C0 2 and form car- 
boxy lie acid, enhancing secondary porosity in pre- 
viously cemented beds by dissolution of calcite 
(e.g:, Schmidt and McDonald, 1979). Hydrody- 
namic forcing by conversion of kerogen to hydro- 
carbons may be on the order of 1CH* s"* 1 (Neuzil, 
1995), or perhaps greater (e.g.. Barker, 1987, 
1990). 

Process Category 3: Mineral Dehydration 
During Diagenesis and Metamorphism 

Dehydration of hydrous minerals during dia- 
genesis and metamorphism is a potential source 
of high fluid pressure (Neuzil, 1995; Hunt, 1996; 
Osborne and S warbrick, 1997). Liberation of 
chemically bound water during alteration of 
smectite to illite is probably not important in 
Cook Inlet basin given the relatively small volume 
of illite in the matrix of lower Tertiary and older 



and dehydration of the crust and subjacent 
emplaced by subduction and underplating 
release fluid to percolate upward beneal 
basin. Cretaceous marine rocks equivale 
those exhumed in the Chugach-Kenai Moui 
accretionary complex were thrust beneati 
plate margin in Late Cretaceous and early Tei 
time (Moore et al. a 1991), and several hun 
kilometers of the Yakutat micropiate was 
ducted during Neogene time (Plafker et a!., 1< 
There is no direct evidence for fluids from 
speculative source, but the broadly distrib 
seismicity at depths between 20 and 35 
(Stephens et al., 1995; Ratchovsky et ah, M 
could be triggered by high-pressure fluid relej 
by dynamo- thermal metamorphism. No 
(1995) estimated T= l(T 15 s" 1 for mineral dc 
dration during metamorphism, a value that 
adopt for comparison with other mechanis 
(Table 3). 

DISCUSSION 

We find that abnormal fluid pressure occi 
sporadically in discrete intervals several meters 
several hundred meters thick in die lower Tenia 
rocks, and speculate that high fluid pressure m 
be even more extensive in older rocks. This lad 
speculation is based on fluid pressure and son 
log measurements in the Swanson River anticlii 
(Fig. 7; Table 1). There is no compelling ev 
dence, however, for a laterally extensive, higl 
. pressure fluid compartment below 4 km ii 
volving Tertiary rocks, nor for a 1 -km-thicl 
calcite-cemented sealing layer that cuts acros 
. .strata and structure at depths between 3 and 4 kn 
as previously proposed (e.g.. Hunt, 1990] 
He-Zhiyong and Banks (1 995) reached a sirnila 
conclusion based on their analysis of fluid pres 
sure in Cook Inlet basin. 

We propose a conceptual model that links Quk 
migration. in Cook Inlet basin to faulting and the 
development of high fluid pressure in Mesozoic 
rocks beneath the basin (Fig. 11). This fault- 
conduit model integrates our diverse geological, 
penological, and geochemicaJ data, together with 
fluid-pressure measurements (Table 1). In this 
model, Na-Ca-Q brine is produced in Mesozoic 
and possibly older basement rocks beneath the 
basin. Fluid is expelled during faulting, migrates 
upward through large reverse- and oblique-slip 
faults, and is injected into the cores of anticlines 
along cross-faults and through folded bedding 
(Fig. 11). Mixing of the Na-Ca-CI brine with 
Na-HC0 3 fluid precipitates carbonate minerals. 
Zeolites are precipitated by Na-Ca-CI fluid reac- 
tion with rock in Mesozoic and lower Tertiary 
strata. Alternating bands of secondary cement re- 
duce permeability normal to bedding and channel 
fluid along permeable interbeds. Fault permea- 
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bility is created by fracturing, and enhanced by 
transport of fluid at high pressure, which in- 
creases fracture aperture by reducing effective 
stress. Conversely, fault permeability decreases 
as fluid pressure dissipates in the underlying 
reservoir (Roberts and Nunn, 1996) and fractures 
become sealed by zeolite and carbonate cement. 
This fault-conduit model is similar to that pro- 
posed by Sibson (1990) for fluid expulsion 
through reverse faults, and by Roberts and Nunn 
(1996) for brine and hydrocarbon migration 
through normal faults. 

The fault-conduit migration model (Fig. 1 1) is 
based on several lines of evidence, which we dis- 
cuss in the following sections. 

Basin Fluids and Secondary Cement 

Processes related to fluid migration must be 
considered in the context of the tectonic setting 
and history of Cook Inlet basin. Sedimentation, 
diagenesis, and tectonics all combined to produce 
a forearc basin in which two chemically distinct 
fluids became vertically stratified because the . 
rate of subsidence was roughly equal to the rate . 
of fluvial sedimentation. Na-Ca-Cl brine that 
originated in marine Mesozdic strata was buried ^ 
beneath Tertiary stream and flood-plain deposits. \ 
Fresh water trapped in the fluvial strata evolved 
into Na-HC0 3 brine during diagenesis. Migra- 
tion and mixing of these two fluids were induced 
by faulting and folding, and ultimately caused 
formation of the zeolite and carbonate cement 
that markedly changed the permeability of the 
lower Tertiary arid Mesozoic rocks. 

The distributions of zeolite and carbonate 
veins and cement are critical aspects of the fault- 
controlled migration model (Fig. 11). Zeolite 
cement is confined to Mesozoic and Paleocene 
strata, but carbonate minerals are found in strata 
ranging from Mesozoic through early Miocene 
age. Zeolities are produced by reaction of Na- 
Ca-Cl brine with Mesozoic sandstone and silt- 
stone. Migration of the Na-Ca-Cl brine into over- 
lying Tertiary rocks proceeds by upward flow of 
fluid along faults, and lateral spreading of the 
fluid away from fault exit points into the Tertiary 
strata (Fig. 11). This migration process forms zeo- 
lite and carbonate mineral veins in faults, and pre- 
cipitates secondary cement in the folded Tertiary 
rocks. The abundant secondary calcite in con- 
glomerate and coarse sandstone beds presumably 
reflects mineral precipitation by transport-limited 
chemical reactions (Fig. 10), where solutes are 
carried to reaction sites by advection (Phillips, 
1991). Initially, the rate of secondary cementing is 
presumably greater in permeable sandstone and 
conglomerate beds where fluid flux is concen- 
trated, and less in interbeds of lower permeability 
siltstone and shale. This process ultimately re- 
duces both the original variation in permeability 
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Figure 8. (A) Permeability versus depth for core plugs in a 30-m-thick interval in the ' 
Formation, Sunfish #3 well, North Cook Inlet anticline. Note that permeability varies i 
with depth and sample grain density. High grain density is caused by carbonate ceme 
plugs pore space and reduces penneability. 1 mD = 1 x 1(H 5 m 2 . 



between adjacent beds, and the mean permeabil- 
ity of the stratigraphic interval (Fig. 8; Phillips, 
1991). These cemented intervals are aquitards to 
fluid flow across bedding; but channel flow 
through permeable interbeds (Fig. 1 1); * 

Geochemical modeling and oxygen isotope 
compositions of calcite provide additional evi- 
dence for the migration of Na-Ca-Cl fluid into 
lower Tertiary rocks. Calcite cement and vein 
minerals in lower Tertiary strata developed both 
by diagenetic reactions between Na-HC0 3 brine 
and rock grains, and by mixing of this brine with 
migrated Na-Ca-Cl brine (Fig. 10). Of 10 calcite 



cement samples collected from lower T 
strata, 6 are relatively light in 6 I8 0 (Fig. 9 
eating that calcite was probably precipita 
reaction between Na-HC0 3 fluid and rocl 
sumably during burial diagenesis. The re 
ing samples of cement have intermedi. 
heavy 5 J8 0, indicating that this calcite pre 
precipitated by mixing of the Na-HC0 3 
and migrated Na-Ca-Cl fluid. This migr 
mode] is also supported by the isotopic co 
sition of calcite veins in faults, most of u 
plot along a 6 I8 0 trend indicative of fluid m 
(Fig. 9). 



TABLE 2. COMPOSITION OF SUBSURFACE WATER FOR GEOCHEMICAL MODELING 


Component 


Tertiary fluid source 


Mesozoic fluid source 




(mg/kg) 


.(mg/kg) 




Evolved meteoric water 


Evolved sea water 


Total dissolved solids 


3755 


19725 


CI 


648 charge balance 


1 1 989 charge balance 


HCO, 


1830 


28 calcite saturation 


Ca 


0.6 calcite saturation 


2809 


Na 


1149 


4603 


K 


2 K-feldspar saturation 


53 K-feldspar saturation 


Fe 


0.04 siderite saturation 


21 siderite saturation 


Mg 


0.01 chlorite saturation 


1 chlorite saturation 


Ai 


1 kaolinite saturation 


0.1 illite saturation 


SiO, 


35 quartz saturation 


31 quartz saturation 
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The importance of faulting in fluid migration 
is evidenced by the widespread occurrence of 
zeolite and carbonate veins in reverse- and 
strike-slip faults (e.g., Fig. 5). The link between 
faulting and fluid flow suggests that those 
processes that create high fluid pressure in the 
Mesozoic basement are most important on a 
basinwide scale. Volumetric strain caused by 
tectonic shortening, and particularly by shear- 
ing and compaction in fault zones, is presum- 
ably most effective in the Mesozoic and older 
rocks (Table 3). Deformation may be aug- 
mented by metamorphism and mineral dehy- 
dration in the middle crust and below, and by 
thermal alteration of organic-rich rock and hy- 
drocarbons. This latter process was more vigor- 
ous in the past than at present, given the hydro- 
carbon maturation history proposed by Magoon 
(1994, and 1998, personal commun.). Sedi- 
mentary loading and compaction in synclines 
may also generate high fluid pressure, perhaps : 
driving fluid laterally into fault zones or upi 
ward along bedding in anticlines (He-Zhiybng 
and Franks, 1995). Thermal alteration of coal 
beds in Upper Cretaceous and lower Tertiary 
strata could also create high fluid pressure, but 
this process cannot be the driving mechanism 
for flow of Na-Ca-Cl brine out of older marine 
strata. 

The fault-controlled fluid-migration model 
provides a rationale for the heterogeneous dis- 
tribution of fluid pressure in some Cook Inlet 
anticlines. The pressure, spatial dimension, and 
lifetime of a plume of injected fluid; depends oh 
the pressure of the source region relative to that 
at the fault exit point, the fluid flux; duration of 
flow, and the permeability and storativity of the 
rocks into which the fluid is injected (Roberts 
and Niinn, 1996). Generation oif long-lasting 
pressure anomalies is facilitated by cata- 
strophic rupturing of the fault conduit, which 
greatly enhances permeability, and therefore - 
fluid flux. The duration of flow is controlled by 
the rate of pressure depletion in the source 
zone, the sensitivity of fracture permeability to 
increases in effective normal stress (Roberts 
and Nunn, 1996), and porosity reduction by 
mineral precipitation. The time scale for cata- 
strophic rupturing is the recurrence interval of 
large earthquakes on active faults, which is 
probably ^lO 3 yr for the large reverse and 
oblique^sip faults coring Cook Inlet anticlines 
(Haeussler et al., in press). Fluid may also seep 
through fault conduits if fracturing occurs by . 
fault creep, as the result of static strain changes 
induced in the basin by great subduction-zone 
earthquakes like the M = 9.2 event of 1964, or 
by dynamic stresses induced as seismic waves 
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Figure 9. Carbon and Oxygen isotopic measurements for calcite in fault veins and cement 
Mesozoic and Tertiary rocks. PDB — Peedee belemnite; SMOW — standard mean ocean wati 
See text for discussion. 

propagate through the basin. Great subducnon- ..». basin (Fig, 2). Some of these folds are known i 
I- zone earthquakes occur at -700 yr intervals < contain abnormal fluid pressure (Table 1 ); othe 
(CombelKck, 1991) Smaller eaMquakes mat^ apparently do not (e g., Middle Ground Shoj 



originate in the subduction zone or lower crust 
beneath Cook Inlet (Bg. I) are felt on a weekly 
* to rhonAly basis. S-./v/'-'V^' '".''r : i 

Migration History and Cenozoic Tectonics 

Fluid migration in Cook Inlet basin occurred in ! 
two primary phases. Hydrocarbons that migrated? 
into traps prior to the onset of Eocene deformation 
were presumably lost during erosion pnor to dep- 
osition of the Hemlock and correlative conglomer- 
ates in late Ohgocene time, according to Magoon 
(1994). Migration of NarCa^ brine out of Meso- 
zoic beds presumably initiated zeolite and carbon- 
ate deposition in Mesozoic and lower Tertiary 
rocks during this early period of deformation and 
fluid iriigratioh. The second phase of oil and brine 
migration followed deposition of the lower and 
middle Miocene Tyonek Formation, and was ac- 
companied by growth of fault-propagation anti- 
clines that form structural traps throughout the 
basin (Figs. 2 and 4; Magoon, 1994). Reverse- and 
oblique-sup faults provided the primary conduits 
for migration of Na-Ca-Cl brine and hydrocarbons 
into the cores of these anticlines, as previously dis- 
cussed. Here we note that siltstone and clay-rich 
layers at the base of coal beds are usually assumed 
to be the hydrocarbon trap rock (Magoon, 1994), 
but in our opinion, secondary cement is probably 
an important, yet overlooked, factor in permea- 
bility reduction and trap formation within Cook 
Inlet anticlines. 

We have studied only a few of the more than 
24 faults and anticlines located within Cook Inlet 



anticline). There are, however, few pressui 
measurements in pre-Tertiary rocks (Table I 
: The fault-controlled fluid migration mode 
(Fig. 11) provides a conceptual framework that i 
amenable to testing and refinement as explo 
ration and drilling continues in the future. Th< 
model has the potential to guide hydrocarboi 
exploration and development strategies, as wel 
as provide insight into possible links betweer 
V tectonic activity, fluid migration, and fluid-pres- 
sure distribution. This latter topic is of consider- 
able importance for understanding fault stability 
and potential for earthquake generation. Careful 
study of the Quaternary geology is required to 
determine which folds are active, and if there is a 
correlation between subsurface fluid pressure 
and active faulting. 

CONCLUSIONS 

Faulting controls fluid migration in Cook In- 
let basin. Fluid is expelled from Mesozoic 
source rocks in the basement, flows upward 
through faults, and is injected into Tertiary strata 
in the cores of anticlines. Fluid migration is pre- 
sumably episodic, with the greatest fluid flux 
following catastrophic fault rupturing and earth- 
quake generation, which may occur every sev- 
eral thousand years in actively growing folds. 
Migration may also occur by fluid seepage 
where fault permeability is maintained by creep 
or stress transients. Fault permeability is re- 
duced or sealed by precipitation of zeolite and 
carbonate veins. 
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kg Mesozoic formation water added to 1 kg Tertiary formation water 




Grams of rock reacted per kg of Tertiary formation water 



Figure 10. Geochernical simulations showing mineral products precipitated during (A) Na-Ca-Cl and Na-HC0 3 brine-mixing reacti» 
(B) Na-Ca-Cl brine-rock reactions; and (C) Na-HC0 3 brine-rock reactions. 
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Hyaroaynamic process 



nyaroaynamic ioruny 
(1/s) 



[Permeability k = 1 £- 18 m 2 ) 



Sedimentary compaction 
Distributed tectonic shortening 
Fault creep and tectonic compaction 
Fautl rupture (Poroelastic) 

Glacial loading 

Hydrocarbon generation and/or alteration 
Prograde metamorphism 



S4 x 10- 15 
-3 x 10r 17 

3x 1(T 14 to9x 10- 13 
Instantaneous pulse of 

1-10 MPa 
1.5 x 10-" to i.5x 10- 12 

*i x i<r 14 

51 x 10- 15 



>5 km 
>100km 
<1 km 

Affects region of several 

kilometers surrounding fault 
<1 km 
52 km 

220 km 



Focused into synclines 
Not reasonable reservoir dimension 
Shearing rate of 1-3 mm/yr 
Pressure pulse decays in periods 

of 10 to 1000 yr 
Pre-Holocene loading 
Thermogenic and biogenic alteration 

of organic matter and hydrocarbons 
Middle crust and deeper source (?) 



Note: Increase or decrease reservior length by a factor of 10 for each decade increase or decrease in permeability. 






Q -Quaternary Sediments 


' Calcite & Siderite Cement Bands 




Tn - Neogene Strata f 


Zeolite & Carbonate Cement 




Tp - Paleogene Strata 


'Na-Ca-Cl Fluid Plume 




Mz • Mesozoic Rock 


High Pressure Fluid Reservoir 




f Neogene Reverse + 
& Transcurrent Faults * 


* Paleogene Faults 



Figure 11. Conceptual model of fluid migration in a fault-propagation fold. High-pressure 
fluid originates in the Mesozoic or older rocks beneath the basin, migrates upward along faults, 
and is injected as plumes of Na-Ca-Cl brine in lower Tertiary rocks. Zeolite cement and veins 
are precipitated by reaction of Na-Ca-Cl brine with sedimentary rock grains. Calcite and other 
carbonates are deposited where Na-Ca-Cl brine reacts during mixing with Na-HC0 3 fluid in 
Tertiarj strata. Some carbonate cement is also precipitated by reaction of Na-HC0 3 brine and 
rock fragments during burial diagenesis. See text for detailed discussion of the migration model. 



The basin contains two different pore fluid 
Na-Ca-Cl brine reacts with sedimentary roc 
fragments to form secondary zeolite cement i 
Mesozoic and Paleocene rocks and veins i 
faults..NarHC0 3 connate fluid in Tertiary strat 
reacts with rock, and mixes with migrate 
Na-Ca-Cl brine to precipitate carbonate mineral 
in sedimentary pore space and fault veins. Flue 
tuations in the partial pressure of C0 2 also resul 
in precipitation of carbpnate cement and veins 
but this process is inefficient relative to fluid-rocl 
and fluid-mixing reactions, which require much 
less volume of fluid per gram of calcite produced 
V y High fluid pressure may originate by a variety 
of mechanisms acting either alone or together at 
different times and locations within the basin. 
Generation of high fluid pressure In the basement 
is presumably important given the evidence for 
fault-controlled' migration of brines and hydro- 
carbons in anticlines! Volumetric strain caused by 
tectonic shortening and shearing, and dynamo- 
thermal metamorphism in the middle crust are 
.. perhaps most important in basement rocks 
Other, alternative mechanisms include alteration 
of organic-rich rocks and hydrocarbons, sedi- 
mentary, loading and compaction, and possibly 
relic pressure generated by glacial loading during 
late Pleistocene time." . ''.[ '" 
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